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An analytical solution for the flow field of a shear flow over a rectangular cavity contain- 
ing a second immiscible fluid is derived. While flow of a single-phase fluid over a cavity 
is a standard case investigated in fluid dynamics, flow over a cavity which is filled with 
a second immiscible fluid, has received little attention. The flow filed inside the cavity is 
considered to define a boundary condition for the outer flow which takes the form of a 
Navier slip condition with locally varying slip length. The slip-length function is deter- 
mined from the related problem of lid-driven cavity flow. Based on the Stokes equations 
and complex analysis it is then possible to derive a closed analytical expression for the 
flow field over the cavity for both the transversal and the longitudinal case. The result 
is a comparatively simple function, which displays the dependence of the flow field on 
the cavity geometry and the medium filling the cavity. The analytically computed flow 
field agrees well with results obtained from a numerical solution of the Navier-Stokes 
equations. The studies presented in this article are of considerable practical relevance, 
for example for the flow over superhydrophobic surfaces. 
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1. Introduction 

Flows in and over cavities are important blueprints for many phenomena in fluid dy- 
namics. On a reduced scale, they illustrate phenomena occurring in various, often more 
complex applications. These might for example be in aer odynamics, where the vortices 
forming inside cavi t ies or grooves in surfaces are studied ( Pan fe Acrivosl 119671: iMoffattl 
19631: iHigdor] 1198.4 IShankar fc Deshnandel lioooh . The correspo nding flow patterns af- 
fect the drag a fluid s treaming over the surface is subjected to (|Shen fc Florvanlll985 ; 
Prasad fc Koselll989h . Traditionally, such troughs, or more generally surface rough- 
ness features, have always been perceived as unwanted, since the vortices forming inside 
the cavities may lead to higher drag or acoustic disturbances. Yet, at higher Reynolds 
numbers, riblets may also reduce drag, which has been of great interest more recently 
( Viswanathl 2002 : Jovanovic et al. 2010h . 

In all these investigations, only a single fluid phase flowing over the cavity and filling it 
at the same time is considered. However, there exist numerous applications, where a fluid 
flows over a cavity filled with a second immiscible fluid. It is well known that by capil- 
lary condensation vapor condenses in sufficiently small cap illaries or channels although 
the vapor pressure is below the bulk saturation pressure (jFisher fc Israelachvililll979t 
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Fisher et g?Jll981r ICiccotti et aLll2008l ). So does vapor inside the roughness features of a 



surface. This may have an immediate impact on the applications in aerodynamics men- 
tioned above, since the condensed water has a viscosity different from the gas and hence 
moves at a slower speed. Moreover, we can also encounter the inverse situation, where 
small bubbles, often of na nometer size, are e nclosed in the roughness features of the wall 
of a water-filled channel (|Yang et aZ.1 120070 . Additionally, al l super hydrophobic, Lotus 
leaf-mimicking surfaces take advantage of such enclo sed air (|Lee et al. 2008; Ma & Hill 
2006; iNosonovskv fc Bhushanll2008t iRothsteinl [2010h . Other examples are the flow over 



membranes or porous media, which may be filled with a second immiscible fluid as well. 

Besides the work on flows over cavities, which includes the fluid motion inside the cav- 
ity, there is also investigation of flows over surfaces with periodically patterned boundary 
conditions. In such a context a cavity is usually represented by a certai n boun dary con- 
dition to the flow above it. The flow inside is not considered. IPhilid (Il972h analyzes 
the flow over a no-shear slot within a no-slip wall in longitudinal and transverse di- 
rection. A no-shear interface is certainly the ideal limit, which cannot be reached by 
a real viscous fluid. Although this condition is often used for modeling air-liquid inter- 
faces (lLauga fc StonefeOO^ISbragagria fc Prosperettil2007t Davis fe Laugall20"09t Crowdvl 
2010t bavies fe LaugalbOig T it is not immediately clear how well air fulfills these ideal 



characteristics. 

The no- shear conditio n corresponds to an infinite slip length 
condition ( Navier 1823h 

I - R du 

U I boundary — P 



in the Navier boundary 



(1.1) 



boundary 



which allows for a nonzero velocity u at the boundary. The slip length j3 relates the 
velocity to its gradient in the direction normal to the surface ^£ an d can be understood 
as an imaginary depth below the surface at which the velocity extrapolates to zero. Still, 
when employing equation (jl.ip for the modeling of a fluid- filled c avity, it is by no means 
clear how to specifically choose /3. Belvaev fe Vinogradova ( 2010l ) have calculated a flow 
over an array of slots with constant slip length. However, it remains to unclear how (3 
is related to the slot geometry or to the properties of the fluid filling the cavities. Fur- 
thermo re, in rea li ty the slip length will not be constant but a spatially varying quantity. 
So far, iHockingj ( 19761 ) has briefly addressed the effect of the viscosity of the medium 
inside a cavity, but only found an approximation to the streamlines in terms of an infinite 
Fourier series in the case of an infinite slip length. 

In this paper, a representation of a rectangular cavity in the form of a boundary 
condition for a second immiscible fluid flowing over a surface is derived. It is directly 
related to the flow pattern inside the cavity, i.e. to the geometry of the cavity and the 
viscosity of the fluid it contains. The effect of these properties on the outer flow can 
therefore be directly explored. 

Based on this, we calculate the flow field of a shear flow directed transversely and lon- 
gitudinally with respect to the cavity. The solution to this mixed boundary value problem 
can be achieved though a representation of the flow by complex variables. Unlike with 
often employed methods based on Fourier series, we can derive a closed-form expression 
and present a simple, explicit formula for the flow field. 



2. Model of the cavity and boundary conditions 

We consider a rectangular cavity of depth h, width b and infinite extension in z- 
direction as depicted in fig. Q] Fluid 1 flows over the cavity, while fluid 2 fills the same. 
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Figure 1: Schematic of the considered cavity with coordinate and velocity directions 



The interface between both fluids is assumed to be flat. For small cavities, this is usually 
justified, since surface tension p rovokes a flat tening of the interface and will only per- 
mit negligibly small curvatures ( Higdonlll985 ). In more general terms it means that we 
consider a flow problem at sufficiently small capillary numbers. 

In the following, a boundary condition is derived, describing the effect of the cavity 
on the outer flow. It can be employed without the need of calculating the detailed flow 
inside the cavity. 

A first-principles modeling of the cavity is achieved by taking into account the vis- 
cosity of the medium inside the cavity and the geometry of the cavity. The modeling is 
exemplified for a flow perpendicular to the cavity. The longitudinal case is analogous. 

At the fluid-fluid interface, continuity of velocity and shear stress applies: 



Ul = u 2 
dui du2 

m-z- = m- 



dy dy 
while at the wall, the no-slip condition holds: 



for \x\<-, y = 0, 



(2.1) 



Ui = for \x\ > 



y = 0. 



(2.2) 



u denotes the velocity in cc-direction and rj is the dynamic viscosity. The indices 1 and 2 
refer to the respective fluids. 

For the velocity u 2 at the upper side of the cavity a Navier slip condition 



U2(x, 0) = 8(x) 



du 2 
dy 



(2.3) 



y=o 



holds, with the local slip length function 8(x) to be determined. Incorporating this equa- 
tion, the two transition conditions (|2.1I) collapse into one boundary condition, such that 
the boundary conditions for fluid 1 are 



Ul = NS(x)^ 
dy 

ui = 



for \x\ < -, y = 
for |ar| > ^, y = 



(2.4) 
(2.5) 
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h/b 0.25 0.50 1.0 2.0 5.0 
d*/b 0.0875 0.1628 0.24 0.25 0.25 



Table 1 : Center of the first vortex after Pan fe Acrivosl ( 1967 ) 



with N abbreviating the viscosity ratio of the two fluids 

N=^-. (2.6) 

m 

Equation (12.41) again has the form of a Navier slip condition. The net slip length for fluid 
1 is thus the product of the viscosity ratio and the slip length function. 

As soon as the slip length function has been determined, the flow over the cavity can 
now be studied without the need of computing the flow field inside the cavity. This is of 
great advantage for analytical calculations. Since in the following the focus is on the flow 
over the surface outside of the cavity, the subscript 1 is no longer needed and is dropped 
from now on. 

A suitable expression for the slip- length function 5(x) needs to be found. At the corners 
x = ±|, obviously 8{x) =0 to make both boundary conditions match. At the center of 
the cavity at x = 0, the flow is least suppressed by the walls of the cavity, and hence, the 
slip length reaches a maximum d. 

Values for d can be obtained from the investigation of li d driven cavitie s. Here , the 



upper boundary of the cavit y is also considered being flat. iPan fc Acrivosl (|1967l ) and 



Shankar fc Deshpande ( 2000t ) have shown that, depending on the aspect ratio of the cav- 



ity, one or more counter-rotating vortices of decreasing intensity develop in the cavity. 
For low Reynolds numbers, the centers of these vortices are always located on the center- 
line of the cavity x — 0. The distance from the center of the first of these vortices to the 
top of t he cavity d* gi ves an upper bound to d. This is due t o the shape of the vel ocity 
field (cf. iHigdonl ( 1985 )) being driven by the upper boundary. IPan fc Acrivosl ( 1967h give 



the center of the first vortex D = d* /b as a function of the aspect ratio h/b (see table [T]), 
which rapidly approaches one fourth of the cavity width as the aspect ratio increases. 

This rather heuristic approach for determining S(x) can be formalized by introducing 
iterative maps 

rj : Ui{x) -> Si(x) (2.7) 

and 

T 2 : Si(x) ^ U i+1 {x). (2.8) 

Pi denotes that from the solution of the generalized lid-driven cavity problem with a 
velocity distribution Ui(x) prescribed at the upper boundary of the cavity a slip- length 
function 8i{x) is determined. T 2 maps a slip- length function used as a boundary con- 
dition at the top of the cavity to the velocity distribution at this boundary. Successive 
application of the map Ti o T 2 generates a series of approximations Si(x), (i = 1,2,...) 
to the slip-length function S(x), i.e. lim^oo Si(x) — S(x), provided that the iteration 
scheme converges. Below we will show that when starting with U\(x) = const., already 
S\ (x) - determined in a simplified manner - yields a good approximation to the flow field 
above the cavity. Therefore, in the framework of this article no further iterations will be 

performed. _ ^ 

In that spirit, we utilize the results of IPan fc Acrivos ( 1967 ) to choose a differentiable 



function S(x) with <5(±|) = and 6(0) = d*. The specific form of 5(x) is still to be 
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determined. This will be done in a way that allows calculating a closed-form solution for 
the flow field, as will become clear in the following section. 



3. Transverse Flow 

In the following, a fluid flow in x-direction over the cavity is investigated and its flow 
field is derived. The fluid is driven by a shear stress Too i n x-direction at y — > oo. 
At sufficiently small Reynolds numbers the flow obeys the Stokes equations 

=0 (3.1) 

=0. (3.2) 

u and v are the velocities in x- and y-direction, respectively, and p is the pressure. With 
the stream function 

-r— = — v and -r— = u (3-3) 

ox ay 

the Stokes equations transform into the biharmonic equation 

AA* = 0. (3.4) 
The solution to th e biharmonic equation is given by Goursat's theorem ( Smirnowlll974l ; 



dp 


d 2 u 


d 2 u 


dx 


dx 2 


dy 2 


dp 


d 2 v 


d 2 v 


dy 


dx 2 


dy 2 



Muskhelishvilil 1 1 9751 ) . a special form of which reads 

* = Re((I- 3 )W t (3)). (3.5) 

The Goursat function Wt(j) is an analytic function to be found, with the subscript t 
indicating the transverse flow direction. Re is the real part and the complex coordinate 
is given by 

3 = x + iy and 3 = x — \y. (3-6) 

In equation (|3.5[) . the general Goursat theorem containing two Goursat functions has 
been reduced to an expression containing only one unknown function by elimination of 
the undetermined constant in the stream function. Due to the impermeabi lity of walls 
lying on the real a xis, the stream function is set to at these boundaries (cf. |Philip1ll972l : 



GarabedianlH9 66). Expressio n (13.51) haspreviously been used in the solution of similar 



problems (|Atanackovidll977[ |Philip|ll972l ). 

With the ansatz (|3.5|) the boundary conditions (|2.4|) and (|2.5|) read 



Im(Wt) = 2N5{x)Re(Wt) for |x| < -, y = (3.7) 

Im(Wt)=0 for N>-, y = 0. (3.8) 

The prime denotes the derivative with respect to 3. 

As an elementary mathematical case, the slip length function 6(x) is assumed to be 
of elliptic form, fulfilling the conditions of section [2j In complex notation 8(x) can be 
written as 

/ b 2 

S(x) = -i2Ddx 2 - — (3.9) 

being real for |x| < 6/2 and imaginary for \x\ > 6/2 (fig. [2]). This form is essential 
for keeping the follow ing calculations an alytically tractable. Nevertheless, one is not 
restricted to an ellipse. Icherepanov ( 19771 p. 68-69) outlined a method of approximating 
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Figure 2: real (solid line) and imaginary (dashed line) part of S(x) 

any arbitrary function being continuous on a finite interval and vanishing at its ends 
according to the same requirements concerning its real and imaginary parts as given 
above. However, for convenievce we will use the expression of eq. (|3.9j) . 

Similarly to the boundary conditions at y = 0, the condition of a given shear stress 
Too = Vi du/dy at infinity can be written in term of W(i) as 

Re(W!) = — for y -> oo. (3.10) 



The solution procedure follows the method of ISotkilava fc Cherepanovl l|l974h . origi- 



nally employed in elasticity theory. An auxiliary function 

Ft® = Wlii) \ ==Wi { i ) (3.11) 

ADN^jf _ £ 

is chosen such that the boundary conditions (|3.7j) and (|3.8j) take the form 

Re(F t ( 3 ))=0 for |ar| < |, y = (3.12) 

Im(F t ( 3 ))-0 for |x| > J, y = 0. (3.13) 



2 



The condition at infinity (|3.10|) is 



Re (^(3)) = ^(l- 5 ^) for y^oo. (3.14) 
The conditions (|3. 121) and p,13p constitute a mixed boundar y value problem for the 



half plane. Its soluti on is given by the Keldysh-Sedov formula (jKeldvsh fc Sedovl 11937 



Muskhelishvilil [20081) . Since (|3.12[) and p,13p are homogeneous, only the homogeneous 



part of the Keldysh-Sedov formula is required, which reads 



In ou r case, L(j) is a linear function with two real coefficients (|Lawrentiew fc Schabat 



19671) to be determined. This is achieved by requiring the condition at infinity (13. 14[) . 



giving 
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The auxiliary differential equation p. lip can now be solved. Its general solution is 

= W (i) (^C + Jf^W^O)^^ (3.17) 

with the solution to the homogeneous problem 



a 



/ 7 -2 d 3 

Waft) = e a ° Wa'-V (3.18) 

and a constant C to be determined. C can be obtained from the condition of symmetry 

lm«( 3 ))=0 at 3 = 0, (3.19) 

ensuring the analycity of W t (z). 

Integrating from 3q = |, i.e. evaluating the Cauchy principal value of the integrals in 
(|3.17l) and (|3.18[) . the result for W t {$) is computed as 



Wt ( 3 ) = ¥f 1. (3.20) 

tW 4771 1 + 4DN V ; 

Following Goursat's theorem (I3.5[) . the expression for the streamfunction is 



^ = 2 ?/ + 4£)ATImW (2; + iu) 2 . (3.21) 



This result is consistent with the work of Philid (1972), who calculated the transverse 



flow over a no-shear slot. In the limit of a very low viscosity fluid in the cavity 



lim * = I» yI m J ( x + iy )2 _ 5! (3.22) 
Af-i-oo 2?7i V 4 

which is equivalent to Philip's result. 

When evaluating ()3.21[) computationally in the second quadrant of z one should be 
aware of crossing a branch cut, which for a square root is usually chosen along the negative 

real axis. To avoid ambiguity one may want to write i\J ^ — 3 2 instead of y^ 2 — ^ for 
a; < 0. 



In fig. [31 the streamlines are plotted in the form 
2r h Y 



* = ^ = 1 + V y + ^W^ + i^) 2 -!] (3.23) 

with X = £c| and Y = y|. The only flow-determining parameters are -D, the ratio of 
the maximum local slip length and the cavity width, and the viscosity ratio N. Fig. [3] 
shows the streamlines for cavities deeper than wide (D = 0.25) and varying TV. Clearly, 
the higher the viscosity ratio, i.e. the lower the viscosity of the cavity medium, the more 
the presence of the cavity influences the flow over the wall. For a high viscosity ratio, the 
flow is notably accelerated towards the cavity, while for a small viscosity ratio the flow 
pattern approaches the one of a plain wall. 

This behavior is explicitly apparent from fig. |4j The velocity u at the boundary y = 
is normalized with Wp !ina x, being the maximum achievable velocity for the perfect slip 
case N^ooatx~y = 0. The velocity obviously increases with N. In the same manner, 
the shear stress at y = is similar to the average shear stress for low N, while being 
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Figure 3: Transverse flow streamline patterns for D = 0.25 and varying viscosity ratio: 
a) N = 100, b) N = 10, c) ./V = 1, d) N = 0.1 
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Figure 4: Normalized a) velocity and b) shear stress along the boundary y = 



reduced towards for high N. It can be further observed, that for the chosen elliptic slip 
length distribution S(x), the shear stress is constant along the cavity. 

It should be noted that for N = 1 the flow field is not identical to the case of a single- 
phase flow. The flow profile is modified due to the flat fluid-fluid interface, which strictly 
separates both fluids and cannot be crossed. 

In the case of water flowing over an air- containing cavity, the maximum velocity at 
x = is Uo/w P ,max = 0.9823 and the shear stress is t/too = 0.0177 for D = 0.25. These 
numbers illustrate the accuracy of previous models, where a vanishing shear stress is 
assumed at the air-water interface. For cavities deeper than wide, the maximum velocity 
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Figure 5: Variation of the velocity at 3 = with the slip parameter D for a range of 
viscosity ratios N 



was hence calculated about 2% too high, which nevertheless, is still a reasonable approx- 
imation. However, for smaller aspect ratios of the cavity, the assumption of vanishing 
shear stress becomes increasingly inappropriate, e.g. for h/b — 0.25, i.e. D = 0.0875 from 
Pan fc Acrivosl (jl967l ). the deviation is about 5%. 

For lower viscosity ratios than for the water-air system (TV ps 56), the influence of 
the cavity depth gets more prominent, which can be observed from fig. [5] Here, the 
dependence of the normalized velocity uo/w Pj max at 3 = on the slip parameter D, which 
is a function of the cavity height, is shown for various values of TV. At N = 1, uq / w Pjma x 
drops down to about 50% when the aspect ratio is lowered from h/b > 1 to h/b — 0.25. 
In the case of air flowing over water (N = 0.018) uo/w Pj max is reduced to approx. one 
third upon the same variation of the aspect ratio. Furthermore, it can be observed, that 
the influence of a variation of N on Mo/^max, translating to a modification of the outer 
flow field, is strongest around N = 1. At very large or small values of N, a variation of 
N has only little effect. This fact is also apparent in fig. [3] (a) and (b), which differ only 
slightly. 

As explained above, S(x) may be a general function with <5(±§) = and a maxi- 
mum 5(0) = d. As a minimum model and in order to keep the calculations analytically 
tractable, we have employed an elliptic form with d taking the value of its upper bound 
known from lid driven cavity flows. In order to assess the accuracy of this model, its 
results are compared to numerical calculations. 

Both the fluid inside the cavity as well as the fluid flowing over the cavity are in- 
cluded in the numerical computation. Exemplary, water (density p = 10 3 kg/m 3 and 
dynamic viscosity r\ = 10~ 3 Pas) and air (p = 1.2 kg/m 3 , 77 = 1.8 • 10~ 5 Pas) have been 
considered. At the interface between the fluids, continuity of velocity and stress (|2.1|) 
holds. The dimensions of the cavity are b = h = 50 urn. At a distance to the wall of 
y = 500 urn, a shear stress of 0.1 Pa is imposed, so that Re = u\ x=0 y=0 bp/rj <C 1 for 
both fluids. The Navier-Stokes equations obeying these conditions have been solved with 
the commercial finite-element software COMSOL Multiphysics (COMSOL Multiphysics 
GmbH, Gottingen, D) on a calculation domain of 3 mm width with a triangular mesh 
of approx. 65000 elements. The mesh was strongly refined along the fluid-fluid interface 
to accomodate the discontinuous change in the boundary condition for the flow over the 
cavity between the fluid-solid and fluid-fluid interfaces. Subsequently, the numerically 
computed velocity was interpolated and integrated to yield the stream function. In fig. 
[6j a comparison between analytical and numerical results is displayed for two different 
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Figure 6: Comparison of the transverse flow streamlines calculated analytically with 
equation (|3.21j) and numerically: a) water flowing over an air-filled cavity b) air flowing 
over a water-filled cavity. Note the different scalings of the y-axis in parts (a) and (b) 



cases: water flowing over an air-filled cavity (a) and air flowing over a water-filled cavity 
(b). In case (a), an excellent agreement between the analytical and the numerical results 
is found. The agreement is still good in case (b), only upon great magnification of a region 
close to the air-water interface some deviations between the analytical and the numerical 
results become visible, c.f. fig|U(b). Specifically, the analytically computed streamlines 
have a higher curvature than the numerically computed ones. This corresponds to a small 
overprediction of the velocity, which is expected owing to the choice of d being an upper 
bound. 



4. Longitudinal Flow 

We now address the longitudinal shear flow over a cavity. The problem is mathemat- 
ically closely related to the transverse flow problem. For an infinitely long cavity, the 
velocity w in z-direction is governed by the Laplace equation 

0. (4.1) 



dx 2 dy 2 

The boundary conditions are equivalent to the transverse case, now applied to the velocity 
component w 

w = NS(x)^ for H<|, y = (4.2) 

w = for |af| > ^, y = 0. (4.3) 

We can seek the solution to the above problem in terms of an analytic function Wi(3). 
The subscript 1 stands for the longitudinal flow direction. By definition, both the real 
and imaginary part of W\{i) fulfill the Laplace equation independently, hence we choose 



w = Im(Wi(3))- 



(4.4) 




Figure 7: Longitudinal flow velocity w/w VtUlsiK for D = 0.25 and varying viscosity ratio: 
a) N = 100, b) N = 10, c) N = 1, d) N = 0.1 

This leads to the boundary conditions 

Im(Wi) = N6{x)Re(W{) for |i| < ^, y = (4.5) 

Im(Wi) = for \x\ > |, y = (4.6) 



2 



and the condition for the stress at infinity 



Re(W{) = — for y -> oo. (4.7) 

Obviously, the structure of this boundary value problem is identical to the transverse 
case (|3.7[) and (|3.8p . Wi(j) directly relates to W(%), since both functions are analytic 
and obey boundary conditions differing only by constant factors. Goursat's theorem can 
therefore also be understood as a relation between the transverse and the longitudinal 
flow. 

Thus, the solution process of the longitudinal case follows exactly the procedure out- 
lined above. This calculation (for intermediate steps see appendix [X| leads to a velocity 



' y + 2DNJm\/ (x + \y) 2 - — I . (4.8) 



7?! 1 + 2LW 

The velocity isolines corresponding to several viscosity ratios and the same cavity geom- 
etry as in the previous section are shown in fig. [71 Analogously, ui Pimax is the maximum 
velocity at y — for perfect slip. We observe the same flow enhancing effect for large 
N as for the transverse case. Still, velocities are higher for t he longitudin al flow, since 
Wp.max = 2itp ima x- This is in accordance with previous results (|Philiplll972h . 
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5. Conclusions 

In this paper, wc have investigated the flow over a rectangular cavity which contains a 
second immiscible fluid. Such a situation may arise in numerous practical applications, for 
example in situations in which a gas flows along along a rough surface with liquid-filled 
roughness features or when water flows along a superhydrophobic surface with entrapped 
gas bubbles. 

In the calculation, the cavity is modeled as a boundary condition to the outer flow. The 
boundary condition appears in the form of a locally varying slip length, the functional 
form of which was determined by referring to previous results from lid-driven cavity flow. 
Up to now, analytical expressions for the flow field over such a cavity were known only 
for some limiting cases like a perfect, viscosity- free fluid filling the cavity or an infinitely 
deep cavity. We have achieved to formulate the problem in such a way that the flow 
over the cavity is directly related to the cavity geometry and the viscosity of the cavity 
fluid. The solution was obtained as an explicit analytical expression for the flow field and 
can be evaluated for arbitrary viscosities and arbitrary aspect rations of the cavity. We 
studied the flow in transverse and longitudinal direction. Both cases are connected via 
Goursat's theorem. The influence of the viscosity of the cavity medium on the outer flow 
was demonstrated, naturally being stronger the lower the viscosity. All results agree well 
with the known limiting cases. The analytical results have been compared to numerical 
calculations. The agreement between both is extremely good. 

We kindly acknowledge support by the Cluster of Excellence 259 'Center of Smart 
Interfaces' and the 'Graduate School of Computational Engineering' at TU Darmstadt, 
both funded by the German Research Foundation (DFG). 



Appendix A. 

The mixed boundary value problem 

Im(Wi) = N5(x)Re(W{) 

Im(Wi) = 
Re(W/) = — 



t i , b 

for \x\ < -, 

for |ar| > ^, 

for y — > oo 



is solved by the method of ISotkilava fc Cherepanovl ([19741 ) 
The auxiliary differential equation 

1 



2DNJf - £ 



transforms the above problem to 
Re(Fi(j)) = 
Im^Gs)) = 
Re(Fi(a)) = — | I - 



y = o 
y = o 



Too 

m 



i 



2DN 



for |a:| < -, 
for |a;| > -, 
for y — > oo 



y = 
y = 



(Al) 

(A 2) 
(A3) 
(A 4) 



The solution for Fi($) is given by the homogeneous part of the Keldysh-Sedov formula 
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(I3.15P . Determining the coefficients according to the boundary conditions yields 

F l (3) = ^(l--^)-=L=. (A5) 
771 2DN I 2 _ 6i 

Equation (IA 1[) can then be solved to give 



Wi ( J ) = Z~f Yi 1_ (A6) 
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